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Abstract 

In this paper, we investigate global dynamics for a system of delay differential equations which describes 
a virus-immune interaction in vivo. The model has two distributed time delays describing time needed for 
infection of cell and virus replication. Our model admits three possible equilibria, an uninfected equilibrium 
and infected equilibrium with or without immune response depending on the basic reproduction number 
for viral infection Rq and for CTL response R\ such that R\ < Rq. It is shown that there always exists one 
equilibrium which is globally asymptotically stable by employing the method of Lyapunov functional. More 
specifically, the uninfected equilibrium is globally asymptotically stable if Rq < 1, an infected equilibrium 
without immune response is globally asymptotically stable if R\ < 1 < Rq and an infected equilibrium with 
immune response is globally asymptotically stable if R\ > 1. The immune activation has a positive role in 
the reduction of the infection cells and the increasing of the uninfected cells if R\ > 1. 

Key words: viral infection; global asymptotic stability; Lyapunov functional; immune response 



1. Introduction 

The mathematical models, based on biological interactions, present a framework which can 
be used to obtain new insights and to interpret experimental data. Many authors have formulated 
mathematical models which describe the dynamics of virus population in vivo and these provide 
advances in our understanding of HIV-1 (human immunodeficiency virus 1) and other viruses, 
such as HBV (hepatitis B virus) and HCV (hepatitis C virus) (see [1-3,5-10, 12, 14-25] and the 
references therein). 

During viral infections, the host immune system reacts with antigen-specific immune response. 
In particular, cytotoxic T lymphocytes (CTLs) play a critical role in antiviral defense by attacking 
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infected cells. To investigate the relation between antiviral immune response and virus load, 
Nowak and Bangham [18] developed the following mathematical model. 



< 




N8y(t)-nv(t), 



qy(t)z(t)-bz(t), 



s - dx(t) - kx(t)v{t) 



kx{t)v{t)-8y{t)-py{t)z{t) 



(1.1) 



where x{t) denotes the concentration of uninfected target cells at time t , y{t) denotes the con- 
centration of infected cells that produce virus at time t , v(t) denotes the concentration of virus at 
time t and z(f ) denotes the abundance of virus-specific CTLs. Uninfected cells are produced at a 
constant rate s and die at rate dx(t). Infected cells are produced from uninfected cells and virus 
at rate kx(t)v(t) and die at rate 8y(t). Free virus is produced from uninfected cells at rate N8y(t), 
where N denotes the total number of virus particles from one cell, and die at rate fXv(t). The rate 
of CTL proliferation is given by qy(t )z(t ) and decay at rate bz{t ) in the absence of stimulation by 
the infected cells. Infected cells are killed by CTLs at rate py(t)z(t). All parameters are positive 
constant. 

Korobeinikov [8] studied global properties of a basic viral infection model which ignores 
immunity ((1.1) with p = 0). By assuming that the incidence rate of infection is given by a 
functional form, more general viral infection models are proposed and investigated (see [6, 9]). 
Wodarz et al. [23] considered a mathematical model for two types immune responses. Murase et 
al. [15] and Kajiwara et al. [7] studied stability of some mathematical models for virus-immune 
interaction dynamics. Recently, Priiss et al., [20] showed that (1.1) always admits an equilibrium 
which is globally asymptotically stable by constructing Lyapunov functions. 

On the other hand, in modeling of many biological processes, time delays are usually intro- 
duced for the purpose of accurate representations of the phenomena. In virus dynamics, it has 
been assumed that new virus particles are produced after the initial infection with a time inter- 
val and this leads mathematical models by delay differential equations. The estimated values of 
kinetic parameters are usually changed by these delay differential equations (see [5, 14, 16, 17] 
and references therein). Mathematical analysis for these models is necessary to obtain an inte- 
grated view for the virus dynamics in vivo. In particular, the global stability of a steady state for 
these models will give us a detailed information and enhances our understanding about the virus 
dynamics. 

In this paper, we introduce distributed (continuous) time delays to (1.1) and study its global 
dynamics. Let h\ and hi be positive constants and/i (t) : [0,/zi] — s> M+ w\Afa{x) : [0,/z2] — > K+ be 
integrable functions with /q 1 fa (z)dT = f£ 2 fa{x)dt = 1. As in Mittler et al., [14] and Nelson et 
al., [17], we assume that the infected cells y(t ) appear after the initial infection with a time period 
T and T is distributed according to /i(t) over the interval [0,/zi], where hi is the limit superior 
of the infection delay. In addition, we assume that a time is needed for the virus production after 
a virions enter a cell (see also [12, 24]). Thus, we also assume the production delay T, which is 
distributed according to fa (t) over the interval [0, hj\, where hi is the limit superior of this delay. 
Then, we obtain the following viral infection model. 
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< 




qy{t)z{t)-bz{t). 



k d / f x (r)x(f - r)v(f - x)dx - 8y{t) - py{t)z{t) 
h 

N d S f 2 f 2 (T)y(t-T)dT-nv(t), 



s — dx(t) — kx(t)v(t), 




(1.2) 



The infection rate k d satisfies k d < k and the total number of virus particles from one cell N d 
satisfies < N, if we incorporate the probability of surviving of the infected cells and virus 
particles between the time for infection and for virus production, respectively. 

Stability analysis for (1.2) with discrete intracellular delay was carried out by Li and Shu [10] 
and Zhu and Zou [25]. Recently, based on Li and Shu [10], Li and Shu [11] has investigated 
a viral infection model with a general target cell dynamics, a nonlinear incidence rate and dis- 
tributed delay. Li and Shu [10, 1 1] showed that their model always admits an equilibrium which 
is globally asymptotically stable and it is necessary to have a logistic mitosis term in the target 
cell dynamics for generating a periodic solution. On the other hand, Zhu and Zou [25] estab- 
lished global stability of an uninfected equilibrium and obtained sufficient conditions for local 
asymptotic stability of two infected equilibria. However, since Li and Shu [10, 11] did not con- 
sider the immune response to the viral infection and Zhu and Zou [25] did not address the global 
stability of the two infected equilibria for their model, the global dynamics of (1.2) is still un- 
clear and, hence, our main aim is to establish the complete global dynamics. We show that (1.2) 
has three possible equilibria, an uninfected equilibrium and infected equilibrium with or without 
immune response and always admits one equilibrium which is always globally asymptotically 
stable. Moreover, it is shown that if the immune response is activated, then the infected equi- 
librium with immune response is globally stable. This implies that the immune response has a 
positive role in the reduction of the infected cells. 

The paper is organized as follows. In Section 2, we show the positivity and ultimately bound- 
edness of the solutions for (1.2) under suitable initial conditions. Then, we introduce two impor- 
tant parameters, the basic reproduction number for viral infection Rq and for CTL response R\, 
defined by (2.4) and (2.5), respectively, and three possible equilibria for (1.2). In Section 3, we 
establish global asymptotic stability of these equilibria by constructing Lyapunov functional. It 
is shown that (1.2) always admits one equilibrium which is globally asymptotically stable and, 
hence, we obtain the complete global dynamics of (1.2). In Section 4, we study discrete delay 
model which was considered by Zhu and Zou [25] and show that two infected equilibria of their 
model is not only locally asymptotically stable but also globally asymptotically stable. In Section 
5, we offer a brief discussion. 

2. Preliminary results 

2.1. Positivity and boundedness of the solutions 

To investigate the dynamics of (1.2), we set a suitable phase space. Let h = max{/ii,/z2}. 
We denote by C = C([— h,0],M) the Banach space of continuous functions mapping the inter- 
val [— h,0] into R equipped with the sup-norm. The nonnegative cone of C is defined as C+ = 
C([—h,0],R + ). From the biological meanings, the initial conditions for (1.2) are 
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x(0) = <pi{0),y(8) = pi(d),v(0) = (P3(e),z(0)=zo for 9 e [-h,0], (2.1) 



where (pi £ C+,i = 1,2,3 and zo > 0. 

Lemma 2.1 Every solution of (1.2) with (2.1) is nonnegative for t > 0. Every solution of (1.2) 
with (2.1) is positive for t > h if zo > and either 

<P2(0) + .fo 1 /l(T)fl>l(-T)^(-T)£/T > 0, or 

'0<P3(0) + .fo 2 /2(T)(P2(-TyT>0. 

Furthermore, every solution is bounded above by some positive constant for sufficiently large 



PROOF. The solution (x(t),y(t),v(t),z(t)) of (1.2) with (2.1) exists and is unique on its maximal 
interval of existence (0, a) for some a > 0. We see that x(t) > for all t 6 (0, a). Indeed, this 
follows from that jjjx(t) = s > for any t E (0,(7) whenx(f) = from the first equation of (1.2). 
It also holds that 

z(t)=zoe I ° iqy{s) - b)ds >0, 

if zo > 0. In particular, z(t ) > if z > 0. 

Let us show the nonnegativity of y(t) and v(f ). Since we have 



(2.2) 



y(t) = ((p2(0)+k d ^ f l (T)x(s-T)v(s-T)dT^(S+pz(u))du ds y-J^ 

v(0 = (^{Q)+N d 8 j'J* 2 f 2 {T)y{s-T)dT^ s ds^e-v, 

from (1.2), y(f) > and v(f) > for t > 0. Now we show y(t) > and v(t) > for t > h, if i) or 
ii) holds. 

First, we assume that i) holds. Suppose that there exists a t\ such that y(t\) = 0. Then, from 
(2.2), 

92(0) +k d [" r f l {T)x{s-T)v(s-T)dT^ {S+pz{u))du ds = Q 1 

Jo Jo 

follows. This leads a contradiction to i). Thus, we obtain 

y(t) >0forf >0. (2.3) 
Next, we suppose that there exists a t 2 > h 2 such that v(f 2 ) = 0. Then, from (2.2), 

q>3(0) +N d 8 f 2 [ k2 f 2 (T)y(s-T)dT^ s ds = 0, 
Jo Jo 



follows. On the other hand, we have 

ft [hi 



JO 



f 2 (x)y(s - x)dxds > 0, for t > h 2 , 



by (2.3). This gives a contradiction. Thus, v(t) > for t > h 2 . Similarly, we see v(t) > for t > 
and y(t) > for t > h\ if ii) holds. 

Now we show the boundedness of each solution. Let 



G(t) = % t h (r)x(t x)dx+y(t) + ? z (t), 
k Jo a 



then we see 
d_ 
dt 



jG(f) =(jJ Q hl h (t) (s - - t) - kx(t - x)v{t - x))d^j 

+ (k d jf *' h (T)x(f - T)v(f - t)^t - 5y(0 - M'Mo) + (MOz(0 - |&(0 

/"*' /! (r)x(f - x)dx - 8y(t) - ?-bz{t). 
Jo q 



sk^ dkd p 
Therefore, it follows that 



jG(t)< S -^-rmn{d,8,b}G(t), 



which implies that x(t),y(t) and z(t) are uniformly bounded on (0,c). Then, v(f) is also uni- 
formly bounded on (0,a). Finally, it follows that (x(f),y(f), v(f),z(f)) exists and is unique and 
positive for any t > h. □ 

Remark 2.2 y(t) andv(t) of {1.2) with (2.1) are identically zero for t > 0, if 

9l{0) = (£3(0) = fl{x)<pl{-x)<f>l{-x)dX = f 2 (T)((>2(-T)dT = 0. 

Jo Jo 

2.2. Possible equilibria 

In this subsection, we show that (1.2) has three possible equilibria. Existence of these equilibria 
is determined by a combination of two threshold parameters 

*o = -V, (2.4) 

a k d N d 

and 

Ri = — „ S , , . (2.5) 

d^— + -^8^ 

°k d N d + k d °q 

Ro and R\ are called the basic reproduction number for viral infection and for CTL response, 
respectively (see Gomez-Acevedo et al. [3]). In particular, Rq denotes the average number of 
secondary virus produced from a single virus for (1.2). 
Theorem 2.3 For ( 1.2), there exist an uninfected equilibrium 

£o = (*o,0,0,0),xo = ^, (2.6) 
an infected equilibrium without immune response 

*-«^.")-(^3(-^)-^(-^)- )- ( "> 

ifRo > 1, and an infected equilibrium with immune response 

p , * * * * A / s b N d 8b 8 (k d N d x\ , 



ifRi > 1. 
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PROOF. First of all, we see that (1.2) always has the uninfected equilibrium Eq. To find other 
equilibria, we consider the following equations 

= s — dx* — kx*v*, 
= k d x*v*-8y*-py*z*, 
N d 8y* uv. 
{0 = qy*z*-bz*. 

Assume that there exists an equilibrium E\ = (x\,y\,v\,Q) withxj > Q,y\ > 0,v* > 0. From the 
third equation of (2.9), we see 

y\ = fig*. (2.10) 

Substituting (2.10) into the second equation of (2.9) gives 
Hence, 

x * = a 

1 k d N d 

Then from the first equation of (2.9), it follows 

s-dx\ _ k d N d ( s d M 



kx\ k\l \ k d N d 

v* and y\ is positive, if Rq > 1 . Consequently, there exists the infected equilibrium E\ if R$ > 1 . 

Next, we assume that there exists an equilibrium Ei = (x|,y|,V2,z|) withxj > 0,y 2 > 0,v 2 > 
().-; 0. We have 

* b „ N d 8 „ 

y 2 = -,v 2 = — — yi, (2.H) 

q jU 

from the forth and third equations of (2.9), respectively. Then, we have 



s s 

H q 



4 -d + kv* 2 - d + k ^ (Z12) 



and 

* _ k d x* 2 v* 2 - 8y* 2 
Z2 — ; > 

PV2 

from the first and second equations of (2.9), respectively. By (2.1 1) and (2.12), we see 
= kjx^N^ _$ = $f k d N d x* 2 _ \ = 5 _ 
Zl P M P A M J P 1 
Thus, is positive if R± > 1 and, hence, there exists the infected equilibrium with immune 
response E-i- Consequently, the proof is complete. □ 

Remark 2.4 For R\ > 1, there exist three equilibria, Eq, E\ andE^. Moreover, we have x\ > x\ 
and y\ > y 2 , since 

4-^ = ^(^-D>o, 

and 

y * y *=*± ( s -d-^—-—-\ = (d-V— + —-\ (R 1 -l)>0 
1 2 kd \ k d N d k d q) k8 \ k d N d k d q) 
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follows. Therefore, for the equilibrium condition, the immune activation has a positive role in the 
increasing of the uninfected cells and the reduction of the infected cells. 

Thus, there exist three possible equilibria depending on the values of Rq and Ry defined by (2.4) 
and (2.5), respectively. We see that Rq > Ry always holds. 



3. Global asymptotic stability of three equilibria 

In this section, we study the global dynamics of (1.2) by employing the method of Lyapunov 
functional. Lyapunov functionals, we construct here, are inspired by McClusky [13] for SIR 
epidemic models with distributed delay. From the following result, we see that (1.2) always 
admits one equilibrium which is globally asymptotically stable and hence, the global dynamics 
of (1.2) is fully determined by Rq and Ry. 

Theorem 3.1 i)IfRQ < 1, then the uninfected equilibrium Eq for (1.2) is globally asymptotically 
stable. 

ii) Assume that either i) or ii) in Lemma 2. 1 holds. IfR y < 1 < Rq, then the infected equilibrium 
without immune response Ey for ( 1.2) is globally asymptotically stable. 

Hi) Assume that zq > and either i) or ii) in Lemma 2.1 holds. If Ry > 1, then the infected 
equilibrium with immune response Ejfor (1.2) is globally asymptotically stable. 
Before giving the proof of Theorem 3.1, we introduce some notations. In the Lyapunov function- 
als, the following function is useful. 



g(x) =x—\ — lnx, forxG (0,+°°). 



g(x) has the global minimum at x = 1 and g(l) = 0. 



For simplicity, we will use the following notation in the proof 




forre [0,h\. 



PROOF, i) We construct the following Lyapunov functional 




(3.1) 



where 



U (t)=k d fi(x) x(s)v(s)dsdx + 8 / 2 (T) / y(s)dsdx. 





We calculate the time derivative of Uo(t) along the solutions of (1.2). We see 
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d 
dt 



x(f) 
xq 



:.V„ ( -1--L ) -C/.VU) tV'/ I \'(/ ) ) 



(dx Q -dx(t)-kx(t)v(t)) 



x{t) 

XQ 



x{t) 



1 f-r I C/X(f ) I — -r-r- 1 



XQ 



x(t) 



xq 



*0 



-fcc(f)v(f) + fcc v(f). 



(3.2) 



Next, we obtain 



dt 



T. b(0 + Trv(0 + 7z(0 



=^ / l /iW-T)v(r-T)dT-5y(0-/jy(0z(0 
Jo 



1 



/I2 



+ — AfcS / / 2 (T)y(* - t)Jt - nv{t) 



■(qy{t)z(t)-bz{t)) 



=k d f 1 MT)x(t-T)v(t-T)dT-8y(t) + 8 [ 2 f 2 (T)y(t-T)dT--^v(t)-^bz(t). (3.3) 
Jo Jo N d q 

Finally, we obtain 

£uo{t) = k d J*' f x (t) (x(t)v(t) - x(t - x)v{t -T))dT + 8 / 2 (t) (y(t) - y(t - t)) dsdx 

= k d (x(t)v(t) - £ /j (T)x(f - T)v(f - t)Jt) + 5 (y(f) - J Q h2 f2(T)y(t - x)d\ 



(3.4) 



Consequently, by adding (3.2), (3.3) and (3.4), we obtain 



dt 



Uo (t) = - k -^x(t)[l- 



k d d 



x(t) 

XQ 



k d d ( xq 



+ k d XQv{t) - £-v(t) - -bz(t) 
N d q 



k d XQ - ) v(f) - -bz(t) 



N d 



x{t) 



= — a —x{t)\ +^-(«b-l)v(0--fe(0<0, for/? <l- 



Hence, every solution of (1.2) tends to Mo, where Mo is the largest invariant subset in j = j 

with respect to (1.2). We show that Mo consists of only the equilibrium Eq. Let (x(t),y(t), v(t),z(t)) 
be the solution with initial function in Mo. Then, from the invariance of Mo, x(t ) = xq and z(t ) =0 
for any t. Now we have ^x(t ) — and hence, it follows v{t ) = for any t, from the first equa- 
tion of (1.2). Then, from the second equation of (1.2), we obtain lim t ^ + „y(t) — 0. Therefore, 
the uninfected equilibrium Eq is globally attractive. Since we have dV ^ < for R Q < 1 and 
Uo(t ) > Uo(t ) — Uo(t), the uninfected equilibrium Eq is stable by Hale and Lunel [4, Section 5, 
Corollary 3.1]. Hence, the uninfected equilibrium Eq is globally asymptotically stable for/?o < 1. 
ii) We construct the following Lyapunov functional 



1 (x(t) 



k d x\v\ 



yjt) 



N d Sy\' 



v(f) 



+ 



k d x\v\q 



z (0 + tfi(0. ( 3 - 5 ) 



where 



JO Jr-r V x l v l / JO J«-T 



dsdx. 



We calculate the time derivative of U\{t) along the positive solutions of (1.2) and show that 

dUj (t) 



dt 



< 0. First, we have 

d 
dt 



x(t) 



x{t) 



(s - dx(t) - kx(t)v(t)) . 



Since s = dx\ + kx\v\ holds, it follows 



d 
dt 



x(t) 



dx(t) 



dx(t) 



dx(i) 



x{t) 
x{t) 

1- 



(dx* + kx\v\ - dx(t) - kx(t)v(t)) 



{dx\ - dx(t ) + kx\ v* - kx(t ) v(t ) ) 

2 

+ fcvl ( 1 



x(f) 



x(t) 



x(t) 



x(t) 
1 



1 



x(t)v(t) 



+ £vt 1 -x f v r - — + v, 



(3.6) 



Secondly, we compute 



d_ 

dt 



yit)\ 
y\ ) 



= f (l - ) fa f Q l fi (T)x(f - x)v(t - x)dx Sy(t) - P y(t)z(t) 
= ^ (i - ) (jf /i W M' - - t) - 5y(0)rfr - M'M')) • 



Since we have 5 = rf I 1 , it follows 



d 
dr 



f) 



= y\ ( l ~ W) ) C h (T) - ™ ? - T) - ) " - ,T (' - J) ) } 



yi 

kdx\v\ 



y(t)J Jo 



x(f - T)v(f - t) y(f ) 



i 



yl 



^Ut-- P y(t)z(t) 



y(t) 



i - 1) jf /t (T) (^-,0^-1 (i - |y) M'MO 



/ 

JO 



-y t + 1 ] rfT - 1 (py(0z(0 - wfc(0) • (3-7) 



yi 



Let us calculate the following 



d 
df 



v(0 

Vl 



^O-l))^! 2 ^-^-^)- 



Since, we have jx — Nd ^7* [ , it follows 



d 
dt 



v(0 



N d 5y 



(3.8) 



Now, we see 

dUijt) 
dt 



phi 

= / Mr) 
Jo 



x{t)v(t) 



11 



*(f - T)v(f - T) 



/ /2(T) 
JO 



y(0 



- / /i(T)(x,V;-ln(x;V r )-x ( , T v r . T + ln(x !jT v ?iT ))c/T+ / / 2 (t) (y, - Iny, -y tT + \ny t . x )d%. 
Jo Jo 

(3.9) 

Consequently, by adding (2.12)-(3.9), we obtain 



dx 



If/ (A = dX ^ 
dt ' k d v\x\ y x{t) 



+Ci(f,T)+C 2 (f), 



(3.10) 



where 



1 



Ci(f,r) = ( 1 -x,v t - j+v t 



+ 



/ /l(T) 
JO 



Xt,xVt,% ■ 



Xt,TV t . x 

yt 



-y t + l)dz + / 2 (t) [y t ,T-v, 



Vt 



+ 1 \dx 



fi (t) (x t v t - In (x f v f ) - x f , T v r , T + In (x fiT v f , T )) c/T 



+ / h{*){yt -Iny, -y ( , T + lny (jT )<iT, 
Jo 



and 



c 2 (0 = -rrL(M0z«-Mz(0)+ 



(3.11) 



fe rf XjV| k d x\v\q \dt 

= -TZfi WW-***)) + ^ ("<*W - f MO) 



(3.12) 



Now we claim C 2 (f ) < for all t > 0. Since we have #i < 1, s < d^ + T d 8 7i holds from { - 2A ">- 
Then 



k8 



s — d- 



10 



q 



(3.13) 



from (2.7) and hence, 

C 2 (f)<0. (3.14) 

For Ci (r,r), it holds that 

+ /i(T)(-ln(^v,)+ln(^, T v fiT ))dT+ / / 2 (t) (-lny, +lny fiT ) 

h 2 



o V *» y» / •>(> V v ( / 



1 1 

l--+ln- 

+ /V) (-^ + 1 +ln^) </r + f-^fi + 1 + dx 

k) V yt yt J Jo \ v, v, j 

- -» (£) - f * ( * fir) rfI - f (f ) dti (3 - 15) 

Consequently, < holds from (3.10), (3.14) and (3.15). Hence, every solution of (1.2) 

tends to M\, where M\ is the largest invariant subset in j dU ^ = oj with respect to (1.2). We 

show that Mi consists of only the equilibrium E\. Let (x(t) ,y(t) ,v(t) ,z(t)) be the solution with 
initial function in M\ , then, it holds that 

*(0 = *!, X(f ~?1 f ~ T) = ^ for almost % e [0,h] and y -^-H = ^1 for almost % G [0,/» 2 ]- 

(3.16) 

From the invariance of Mi, we have Jrx(0 = and it then follows that v(f) = v* for any t 
from the first equation of (1.2). From (3.16), we obtain y(t) — y(t — t) = y\ for any t and then, 
z(t) = follows from the second equation of (1.2). Therefore, the infected equilibrium without 
immune response E\ is globally attractive. Since we have dU l^ < and U\ (f ) > Ui(t) — U\ (t), 
the infected equilibrium without immune response E\ is stable by Hale and Lunel [4, Section 5, 
Corollary 3.1]. Hence, the infected equilibrium without immune response E\ is globally asymp- 
totically stable for/?! < 1 < R Q . 

iii) We construct the following Lyapunov functional 

w " k (f ) + 4*' (f ) + (f ) + ^ (f ) + ^ 

(3.17) 

where 

Jo Jt-T \ x 2 v 2 J Jo Jt-t \ y-i j 

Similar to (2.12), we obtain 

f x(t) 



d 
dt 
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We also obtain 



d_ 

dt 



Since we have 8y\ = k d x* 2 v* 2 — py* 2 z 2 , it holds 



S = —(k d x* 2 v* 2 -py* 2 zl). 

12 



Then 

d 
dt 



y\ J 



i 



y(t) 
yl V y(t) 



1 - - ] I k d jf ' f x (T)x(f - T)v(t - T)dT - 1 (^vj - W2Z2)y(0 ~ M'MO 



*'/i(t) f^x(f - x)v(t - T)-k d x* 2 vf-^ ) <lr t(,K 2 y(D - pyU ):.(/ ) ) 
o \ y 2 



>2 

1 



K 1 - 1) ) r 7i (t) <y - T)v(f - t) - ) ^ 



y(0 



(pz2y(0-w(0z(0) 



— -I)(f-1f 



= (l -£) ^h{x){^- Vt ^y t )dx +P zl (l - I) fo-y*) 



+ 1 </T+/>z5&-i)(i-a). 



Similar to (3.8), we also obtain 



d 
dt 



v(0 



(3.19) 



(3.20) 



Let us calculate 



z(0 
z| 



1 (l - Jy) («y(0z(0 - MO) = ^ (i - Jy) («y(0z(0 -©fcW) • 

^-iW)-^) 

Z 2 / 



= ^2 (Zr-1)(>' ( -1). 

Similar to (3.9), we obtain 



(3.21) 
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dU 2 {t) 



rh _ _ f h 2 

= Jo /i( T )( x ' v '~ ln ( x ' v ') _x '>TV ( ,T + ln(x r:T v r:T ))d/T + J f 2 {x)(y t ~\ny t -y tz + \ny tx )dx. 



Consequently, by adding (3.18)-(3.22), we obtain 



(3.22) 



(3.23) 



where 

C 3 (f,T) 



1 -x t v t - = +v t 
x, 



+ jf*' f x (t) {x t ^ x -y t - + l^jdx + £ 2 f 2 (x) (y Xx -v t - } ^ + ^jdx 



+ 



y -hp£ iy, i) (i - *) + t%t<b\ % - 1) (y* - 1) 



k d x* 2 v\ k d x* 2 v* 2 q ' 

f h[ __ _ f hl 

+ / fi(x)(x t v, -\n(x t v t ) -x t:t v,^+ln(x t:t v t , x ))dx + / f 2 (x) (y t -\ny t - y t T +lny t t ) dx 
Jo Jo 



1 -x ( v, - = +v t 
x, 



+ jf*' f x (t) (x t , T v t , T -y t - + l^jdT + £ 2 / 2 (t) (y f , T - v, - ^ + 1) Jt 

f hl _ __ _ f hl 

+ / /i(T)(x,V(-ln(x,v,)-X(, T v (jT + ln(x(, T v r , T ))c/T+ / / 2 (t) (y, -lny, -y, T + lny ; T ) c/t. 
Jo Jo 



Similar to (3.1 1), we see 



-CfM^) dx -C ms (!f) dx ^ (3 - 24) 

Thus, < holds from (3.23) and (3.24). Hence, the solution of system (1.2) limit to M 2 , 

where M 2 is the largest invariant subset in j dU ^ = j with respect to (1.2). We show that M 2 

consists of only the equilibrium E 2 . Let (x(t),y(t), v(t),z(t )) be the solution with initial function 
in M 2 , then it holds that 

X (t) = x * : x(f ~ T | V | f ~ T) = ^ for almost x € [0,hi] and = ^ for almost x e [0,A 2 ]. 

X 2 V 2 ?2 ^2 V 2 

(3.25) 

From the invariance of M 2 , we have ^x(t ) =0 and it then follows that v(f ) = for any f from 
the first equation of (1.2). From (3.25), we obtainy(f) =y{t — x) ^y^foranyf and then, z{t) ~z 2 
follows from the second equation of (1.2). Therefore, the infected equilibrium with immune re- 
sponse E 2 is globally attractive. Since we have dU j^ < and U 2 (t) > U 2 {t) — U 2 (t), the infected 
equilibrium with immune response E 2 is stable by Hale and Lunel [4, Section 5, Corollary 3.1]. 
Hence, the infected equilibrium with immune response E 2 is globally asymptotically stable for 
fli > 1. 
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Finally, the proof of this theorem is complete. □ 



4. Applications 

Our approach is applicable for discrete delay models. Zhu and Zou [25] studied the following 
viral infection model with cell mediated immunity. 
- d 

— x(t) = s — dx(t) — kx(t)v(t), 
a 

—y(t) = ke- Sr x(t - t)v(* - t) - 8y(t) -py(t)z(t), 
d J (4.1) 
—v{t)=N8y(t)-nv(t), 

-z(t) = qy(t)z(t)-bz(t), 

with the initial conditions x(9) = <pi(0),y(O) = y Q ,v(9) = <p 3 (0),z(O) = zq for 9 e [-T,0], 
where (pi(9) £ C([— T,0],K.+),/ = l,3,yo > and zq > 0. All parameters are positive constant. 

For (4.1), similar to (1.2), there exist three possible equilibria. From (2.4) and (2.5), the basic 
reproduction number for viral infection and for CTL response are given by 

$ s 

R = tj — and R\ = 



respectively. There exist the uninfected equilibrium Eq(xq,0 ,0 ,0) ,xq = 4, the infected equilib- 
rium without immune response E\ (x\y j , V\ , 0) (x\ , y i , v\ > 0) if Rq > 1 and the infected equilib- 
rium with immune response E2(x2} '2^2, Zi) (x2,y2^2,Z2 > 0) if R\ > 1 (see also [25, Section 
3]). 

Zhu and Zou [25] established the global asymptotic stability of the uninfected equilibrium £b 
for Rq < 1 . Moreover, they obtained sufficient conditions for the local asymptotic stability of 
infected equilibria E\ and E2 by analysis of associated characteristic equations. Complete global 
dynamics for (4.1) is not clear and an open problem. However, similar to Theorem 3.1 in Section 
3, we establish the following result. 

Theorem 4.1 i) I/Rq < 1, then the uninfected equilibrium Eq for (4.1) is globally asymptotically 
stable. 

ii) Assume y Q + /q 1 /i(t)(Pi(-t)(^(-t)Jt > 0, or (p 3 (0) > 0. IfR\ < 1 <^o> then the infected 
equilibrium without immune response E\ for (4.1) is globally asymptotically stable. 

Hi) Assume zo > and either yo + J /i '/i(t)^i(-t)93(-t)c/t > 0, or ^3(0) > 0. If~R~i > 1, 
then the infected equilibrium with immune response E2 for (4.1) is globally asymptotically stable. 
Zhu and Zou [25, Theorems 3.3, 3.4] showed that the infected equilibrium without immune re- 
sponse E\ is locally asymptotically stable for R\ < 1 < Rq and the infected equilibrium with 
immune response E 2 is locally asymptotically stable for Ri > 1 if the intracellular delay T sat- 
isfies a condition (see [25, Theorem 3.4]). However, by Theorem 4.1, we establish that E\ is 
not only locally asymptotically stable but also globally asymptotically stable for R\ < 1 < Rq. 
Moreover, E2 is globally asymptotically stable, whenever it exists, that is, R\ > 1. 

5. Discussion 

In this paper, we study global dynamics of delay differential equations for a virus-immune 
interaction in vivo. Two distributed time delays represent the time needed for infection of cell 
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and virus replication. Stability analysis for (1.2) with discrete intracellular delay was carried out 
by Li and Shu [10] and Zhu and Zou [25]. Li and Shu [10] studied a viral infection model which 
ignores the immune response to the viral infection and showed that their model always admits an 
equilibrium which is globally asymptotically stable. Recently, Li and Shu [11] has investigated a 
general viral infection model with distributed delay which also does not incorporate the immune 
response. Zhu and Zou [25] established global stability of an uninfected equilibrium and obtained 
sufficient conditions for local asymptotic stability of two infected equilibria when the distributed 
delay in (1.2) is given by a discrete. Zhu and Zou [25] did not address the global stability of two 
infected equilibria for their model. 

To obtain an integrated view for the virus-immune interaction dynamics in vivo, we investi- 
gate the global stability of (1.2) by employing the method of Lyapunov functionals which are 
motivated by McClusky [13] for delayed epidemic models. (1.2) has three possible equilibria, an 
uninfected equilibrium and two infected equilibria with or without immune response. A combi- 
nation of the basic reproduction number for viral infection Rq and for CTL response R\, defined 
by (2.4) and (2.5), respectively, determine the existence of these equilibria. Moreover, they also 
fully determine the global dynamics of the model. The uninfected equilibrium Eo is globally 
asymptotically stable if Rq < 1 and the viruses are cleared. The infected equilibrium without 
immune response E\ is globally asymptotically stable if R\ < 1 < Rq and the infection becomes 
chronic but with no persistent immune response. The infected equilibrium with immune response 
E2 is globally asymptotically stable if R\ > 1 and the infection becomes chronic with immune 
response. Theorem 3. 1 is an extension result of the global stability results in Priiss et al., [20] and 
Li and Shu [10]. Moreover, we improve stability results in Zhu and Zou [25] (see Section 4). 

We see that virus eventually persists if Rq > 1, because the infected equilibrium E\ or E2 is 
globally asymptotically stable in this case. The infected equilibrium without immune response 
Ei is globally asymptotically stable and the immune response does not work for Ro > 1 > R\. 
On the other hand, the immune response is activated and there exist two infected equilibria E\ = 
(x\ , y\ , v j , 0) and E2 = {x\ , y\ , v\ , z\ ) for R 1 > 1 . Moreover, in these equilibria, one can see that the 
relations x\ < x\ and y* 2 < y\ hold due to the effect of immunity (see also Remark 2.4). Therefore, 
the global stability of E2 f° r R\ > 1 indicates that the immune activation has a positive role in 
the reduction of the infected cells and the increasing of the uninfected cells for/?! > 1. 
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